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Abstract 

A Lienard type nonlinear oscillator of the form x + kxx + ^-x 3 + X±x = 0, which may also be 
considered as a generalized Emden type equation, is shown to possess unusual nonlinear dynamical 
properties. It is shown to admit explicit nonisolated periodic orbits of conservative Hamiltonian 
type for Ai > 0. These periodic orbits exhibit the unexpected property that the frequency of 
oscillations is completely independent of amplitude and continues to remain as that of the linear 
harmonic oscillator. This is completely contrary to the standard characteristic property of nonlinear 
oscillators. Interestingly, the system though appears deceptively a dissipative type for Ai < does 
admit a conserved Hamiltonian description, where the characteristic decay time is also independent 
of the amplitude. The results also show that the criterion for conservative Hamiltonian system in 
terms of divergence of flow function needs to be generalized. 

PACS numbers: 05.45.-a, 02.30.1k, 45.20.Jj, 45.05. +x, 82.40.Bj 
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I. INTRODUCTION 



Nonlinear oscillator systems are ubiquitous and they model numerous physical phenomena 
ranging from atmospheric physics, condensed matter, nonlinear optics to electronics, plasma 
physics, biophysics, evolutionary biology, etc. |H Q El U, 0]- One of the most important 
characteristics of nonlinear oscillations is the amplitude or initial condition dependence of 
frequency for nonisolated periodic orbits . For example, for the cubic anharmonic 

oscillator 

x + lu 2 x + (3x 3 = 0, (cu 2 Q ,p>0), (1) 

where dot denotes differentiation with respect to t, the general (periodic) solution is 
x(t) = A cn(Qt + 5), where A and 5 are arbitrary constants with the modulus squared 
of the Jacobian elliptic function m 2 = 2 (S+pA^) anc ^ frequency ^ = a/wq + ft A 2 . In the case 
of limit cycle oscillations the dependence of the initial condition is manifested in the form 
of suitable transient time to reach the asymptotic state. For chaotic oscillations, of course, 
there is the sensitive dependence on initial conditions. In this paper we identify a physically 
interesting and simple nonlinear oscillator of the Lienard type, which is also a generalized 
Emden type equation, that admits for a particular sign of the control parameter nonisolated 
conservative periodic oscillations, exhibiting the remarkable fact that the frequency of os- 
cillation is completely independent of amplitude and remains the same as that of the linear 
harmonic oscillator, thereby showing that the amplitude dependence of frequency is not 
necessarily a fundamental property of nonlinear dynamical phenomena. We also show that 
rewriting the underlying equation of motion as a system of two first order coupled nonlinear 
differential equations the basic criterion for conservative system in terms of the divergence 
of the flow function has to be generalized. 

We consider a nonlinear oscillator system of the form 

k 2 

x + kxx + —x 3 + \ix = 0, (2) 
9 

which is of the Lienard type 'x + f(x)x + g{x) = 0, where f(x) = kx and g{x) = ^-x 3 + Aix 
in the present case. Eq. (J2J) can also be interpreted as the cubic anharmonic oscillator 
defined by Eq. (JTJ) (with u 2 = Ai and (3 = but acted upon on by a strong damping 
type nonlinear force kxx. This type of equation is also well studied in the literature for 
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almost two decades as a generalized form of Emden type equation occurring in the study 
of equilibrium configurations of a spherical cloud acting under the mutual attraction of its 
molecules and subject to the laws of thermodynamics Eq. (J2J) is a special case of a 
eneral second order nonlinear differential equation possessing eight Lie point symmetries 
3] and that it is linearizable. In particular the special case Ai = has been well analyzed 
by Mahomed and Leach []| and explicit forms of the eight symmetry generators satisfying 
sl(3,R) algebra have been well documented. The latter case (Ai = 0) has the exact general 
solution 



where I\ and 1% are the two integrals of motion with the explicit forms 

'1 = -'+!^ (4a) 

h = -/ + (4b) 
o |x 2 + X 

Obviously for Ai = case the initial value problem of Eq. (j2J) appears to be a dissipative 
nonlinear system. 

In the following we show that the nonlinear oscillator equation (J2J) possesses several 
unusual and interesting features, which are contrary to the standard characteristic features 
associated with usual nonlinear oscillators. For example, in the case Ai > we prove that 
the system (j2J) possesses oscillatory motion whose frequency is completely independent of 
the amplitude. Thus we prove that the standard notion of amplitude dependent frequency of 
oscillations is not a necessary condition for nonlinear oscillators. Moreover, we also point 
out that the system (J2J) admits a Lagrangian formulation and from which a conservative 
Hamiltonian can be obtained. Further, we also find that the standard test (divergence of 
the flow function) to identify whether the given system is conservative or not fails here. 
To overcome this situation we propose that the necessary condition for a system to be 
a conserved Hamiltonian is that the time average of the divergence of the flow function 
should vanish instead of the flow function itself vanishing. More interestingly, we also 
find that the system (J2J) when Ai < also admits a Lagrangian as well as a conserved 
Hamiltonian description similar to the case Ai > 0, eventhough the general solution shows a 
dissipative/damping/frontlike aperiodic structure, depending on the initial condition. The 



divergence of the flow function for the case Ai < can be negative for all times in spite 
the conserved Hamiltonian nature. Again to redeem the situation, one requires the time 
average of the divergence of the flow function to be vanishing. Finally, we also point out 
that the special features of the system can be traced to the existence of certain linearizing 
and canonical and nonlocal transformations, which points out to generalization of the results 
in different directions. 

The plan of the paper is as follows. In the following section, we utilize the so called 
modified Prelle-Singer procedure and derive two functionally independent integrals of motion 
for the Eq. (j2J. In Sec. we show that for the case Ai > the system (J2J) admits periodic 
solutions and construct the associated Lagrangian as well as Hamiltonian for this system. 
We discuss the unusual features associated with the system in Sec. |TV] and point out the 
amplitude independence of the frequency of oscillations as well as the nonzero value of the 
divergence of the flow function in spite of its conservative nature. We analyse the nature 
of solutions in the regimes Ai = and Ai < in Sec. El and point out that in spite of the 
decaying/damping nature or frontlike aperiodic nature of the solutions, the system admits 
Lagrangian and Hamiltonian descriptions. In Sec. IVI[ we linearize the Eq. (J2J) through 
different kinds of transformations (invertible point, nonlocal and canonical transformation) 
to trace the unusual features of the system. In Sec. IVIH we generalize Eq. (J2J and compare 
the dynamics with certain other interesting nonlinear oscillators. Finally, we present our 
conclusions in Sec. IVIIII 



II. EXTENDED PRELLE-SINGER PROCEDURE 



Now let us consider the general case Ai 7^ 0. One can proceed to solve the equation 
explicitly by the so called modified Prelle-Singer (P-S) method 0, 0] which identifies the 
integrals of motion and explicit solution, if they exist. One can also use other methods as 
well; however we find that the P-S method is quite convenient |9] to obtain both the integrals 
of motion and solution explicitly. Assuming the existence of an integral I = I(t, x, x) for 
Eq. (J2J) and rewriting the latter as x = P(x,x), where P = — kxx — ^-x 3 — \\X, so that 
Pdt — dx = 0, and introducing a null term S(t, x, x)xdt — S(t, x, x)dx, we find that on the 
solutions, (P + Sx)dt — Sdx — dx = 0. Since I is a constant of motion, it follows that 

dl = I t dt + I x dx + I±dx = 0, (5) 
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so that 



dl = R(P + Sx)dt - RSdx - Rdx = 0, 



(6) 



where R is an integrating factor. Comparing equations (jSJ) with © we have, on the solutions, 
the relations 

I t = R(P + ±S) 
I x = -RS 

I± = —R. 



Then the compatibility conditions, I tx = I x t, h± — Ixt, Ixx — I±x, require that 

D[S] = -P X + SP X + S 2 , 
D[R] = -R(S + P X ) J 



(?) 



where 

D — ® + x ® + P ® 
dt dx dx 

Solving (JZj) systematically for S and R, one can write down the form of the integral of 
motion, 



I = n — r 2 



R+ Tx {n - r2) 



(8) 



where 



J R(P + xS)dt, r 2 = J (RS + —rjdx, 



for the given form P = —kxx — ^-x 3 — X\x. We find two independent sets of compatible 
solutions of Eq. (JJJ), namely, 



—x + |x 2 
S\ = , R\ 



xe 



-2-J-Xit 



X 



(x + \x 2 



S- 2 



2 k 



ix 



\x + V-Ai 



R2 



(fx + V ZI ^)e V ^ XIt 

(x + |x 2 + y/-\ix) 2 ' 



(9a) 
(9b) 
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Consequently, substituting R and S in (|SJ), we find the two (time dependent) integrals of 
motions for Eq. ^ as 

^4^( ^f + ft\ dob) 



III. PERIODIC SOLUTIONS, LAGRANGIAN AND HAMILTONIAN DESCRIP- 
TION FOR THE CASE Ai > 

As we mentioned in the introduction Eq. (J2J) admits two different kinds of dynamics 
depending on the sign of the parameter Ai. In this section we discuss the case Ai > 0. 

A. Periodic solutions 

Now we note that for Ai > both the integrals I\ and I2 become complex. To identify 
two real integrals, we consider the combinations 

4 (3i + kx 2 ) 2 + 9u 2 x 2 
Jl ~ kHJ 2 ~ {3kx + k 2 x 2 + 9u 2 ) 2 ' ( ' 

where u> = a/Ai and 

2e iS i(wt+8) ( 3x + kx 2 - Mux \ 
2 ' k\hl 2 \ \3kx + k 2 x 2 + 9u 2 J ' 1 1 

so that 3i and 1 2T 2 1 can be taken as the two real integrals of Eq. (J2J) for Ai > 0. In Eq. f)12|) . 
5 is a phase constant. 

Thus for the case Ai > 0, the integrals (fTTj) and (JT2J) lead to the explicit sinusoidal 
periodic solution (Fig. 1), 

A sin (cut + 5) 
.-(^Acosiut + Sy 

Q 

< A < — , u= y^, (13) 



x(t) 



where A = 3u\/3i and S is an arbitrary constant. The form of the first integral (jllj) also 
establishes that the system is of conservative type. Note that in (JT3j) 

A 



and 



so that 
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A 



B 



Hat 

the solution ()1HJ) becomes singular at finite times. 



and B may be called as the amplitude of oscillations in the present case. Also for A > 
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FIG. 1: (Color online) Harmonic periodic solutions of the nonlinear Lienard type oscillator for 
different initial conditions (Eq. ((21)) with frequency same as that of the linear harmonic oscillator. 



B. Lagrangian and Hamiltonian description 

The form of the time independent first integral 3i as given by Eq. suggests that one 
can indeed give a Lagrangian and so a Hamiltonian description for Eq. (0) when Ai > 0. 
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Indeed we identify a suitable Lagrangian for the system (J2J as 



2P Ux + f x 2 + 3A J 2A; X 2P ' U4J 



3 

Then the canonically conjugate momentum can be defined as 

= dL 27Xf ( 1 V 3A X 

P 9x 2A; V(A:i: + fx2 + 3A 1 )V 2A;' 1 j 

Consequently, one can obtain the Hamiltonian associated with Eq. (J2J) as 

2 V(A;x + f x 2 + 3Ai) 2 / 2 ' { } 

or 

9A 2 / 2kpi k?x 2 2kp 2fcVp \ 

H " " 2(1 " 3a7 } + ^a7"3aT~^taT ) - ( } 

Note that in the above expressions L,p and iJ, we have kept certain prefactors and constant 
terms in order to have the correct limits for k — > 0, namely the harmonic oscillator limit. The 
nature of the trajectories in the (x — x) and (x — p) planes are shown in Figs. 2. They form 
closed concentric curves in the region |x| < %. 
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Finally we may point out that because of the form of the second and third terms in the 
right hand side of $H]) one may define a modified Lagrangian 

1 



kx + f x 2 + 3Ai ' 

(17a) 



with the conjugate momentum 

k 



P 

//.•I--!- 



{kx + fx 2 + 3A X ; 



(17b) 



and the associated Hamiltonian 

2kx + ^x 2 + 3Ai 



H 



{kx + ^a; 2 + 3A^ 2 



(17c) 



or 



(17d) 

which again correspond to the equation of motion (J2J), provided k ^ 0. However in the 
limit & — > 0, in order to get the correct harmonic oscillator it is preferable to have the 
forms (|T^|) - (fTBj) . On the other hand, the later forms (fTTj) have the advantage that they give 
interesting Ai = limit, which we discuss in the next section. Consequently, we have chosen 
to present both the forms here. 



IV. UNUSUAL FEATURES OF EQ.© WITH Ai > 
A. Amplitude independent frequency of oscillations 

The remarkable features of the general solution (|T3*j) are that the nonlinear system (J2J) 
admits harmonic periodic solutions (|13|) and that the frequency of these periodic oscillations, 
v = y~ = is completely independent of the amplitude or initial condition, a feature 
which is generally quite uncommon for nonlinear oscillators of conservative type. The form 
of the solution is shown in Fig. 1. It may be noted that in the limiting case of vanishing 
nonlinearity, k = 0, one recovers the solution of the linear harmonic oscillator as it should 
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be from Eq. (fT3*j) and there is no change in the angular frequency u = \/Ai even when k ^ 0, 
as can be seen from Eq. However, it may be noted that in the limit Ai or cu —>■ 0, 

A — > and so x — > 0, the trivial solution to (j2J and not to the aperiodic bounded solution 
(JHJ). Thus at Ai = 0, a bifurcation occurs. 



B. The nature of flow function 



Another interesting feature of Eq. (J2J) is that when written as a system of two first order 
equations, it takes the form 



x = y = h{x,y), 
k 2 

y = -kxy - —x 3 - Xxx = f 2 (x, y). 



(18a) 
(18b) 



For Ai > 0, it has one equilibrium point (0,0) which is of centre type compatible with 
the occurrence of nonisolated periodic orbits. Interestingly the conventional criterion for 
conservative systems (see for example Refs. 11, Q 13 ] ) , namely, the value of the 

divergence of the flow function A = ^ + ^ for Eq. (|T8|) vanish, fails here as A = — kx. 
Only the average of A, namely 



A 



1 

T 



Adt = 

2vr 



on actual evaluation using the solution (|13|). vanishes. Such a generalized criterion that 
the average of the flow function vanishes for a conservative system seems to be a necessary 
condition as the present example shows. 



V. THE APERIODIC CASES Ai < 

In the previous section, we discussed the dynamics of (J2J for the case Ai > 0. In the 
following we explore the dynamics of the system (J2J with Ai < 0. 

For Ai < 0, both the time-dependent integrals ()10j) are real, from which the explicit 
solution can be obtained straightforwardly as 

/ 3VW(/i^-l) \ (19) 
VA;/ 1 / 2 eV^' + fc(l + / 1 e 2 v^*)/ 
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where l\ and I2 are constants. Solution (fT§j) clearly shows the dissipative/damping/aperiodic 
nature of the system for Ai < in Eq. (j2J). 

We now note that system Ai < still admits the Lagrangian and Hamiltonian descrip- 
tions. In fact the Lagrangian L and the Hamiltonian H for the case Ai > 0, namely the 
forms (|14|) and respectively, are valid here also. The only difference here is that one 
has to replace Ai with — |Ai| in the respective places in Eqs. dHJ) and (fTBj) for the present 
case. On the other hand the forms (j!7|) are valid for Ai > as well as Ai < 0, though not 
valid in the limiting harmonic oscillator case k — > 0. 

The major difference in the dynamics comes from the nature of the solution. In the present 
case it admits decaying type or aperiodic (front like) solutions only, see Fig. 3, where the 
time of decay or approach to asymptotic value is independent of the amplitude/initial value, 
which is once again an unusual feature for a nonlinear dynamical system. 




2 4 6 8 10 12 14 

t 

FIG. 3: (Color online) Solution plot for case Ai < for different initial conditions. 

The nontrivial general solution for the limiting case Ai = is given in Eq. ©. Interest- 
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ingly in this case also one can deduce a Lagrangian and the associated Hamiltonian, which 
(from eqs. (|I7j)) turns out to be 

L=( - rg A (20) 

and 



/■j ^y> _l_ ,y>2 / 



fc 2 o /-P 



H=-*>,-2j-£, (21) 



where 

P ~ ~(H; + fx 2 ) 2 ' ^ 
The above form of Lagrangian can also be deduced from the form (I14j) . by taking the 
Lagrangian for the limiting case Ai = as 



.. 2k 2 3Al 9A 2 , 

hm — 7^{L —x-\ — 

Ai-o 27Af 2k 2k 2 ' 



The expressions p and H given above can be obtained from this L. The solution plot 
(corresponding to solution (jSJ)) is given in Fig. 4. 

One can also calculate the flow function in the case Ai < 0. Again from Eq. (J18)) we 
find that the divergence of the flow function for Ai < is A = — kx. But from the general 
expression for the solution, one can always choose x(t) > for all finite t. So one might 
conclude that the system is dissipative (because A < 0) inspite of its just proved conservative 
nature. In order to overcome this incompatibility, and noting that for the case Ai < 0, the 



dynamical variable x(t) asymptotically tends to a constant value x(oo) = jt^-, we again 
define the condition for conservative flow is that the time average of the divergence of the 
flow function with reference to the asymptotic (t — > oo) value should vanish: 

A = lim^ i r o T [A(t) - A(oo)]dt = 0. (23) 

In the present cases, we have 

A = - limr^oo | £[x(t) - x(oo)]dt. (24) 

For the case Ai = 0, from the solution (JHJ) we note that x(oo) = and f Q T A(t)dt = —3 log(l + 

fcT( 6/2 +T) ) so that ^ = from ®- 0n the other hand for Al < °> So A (t)dt = -3^f\{T - 
3 log( Il+Ill2 ^ 1+Il+ 1 I ^ -) and A(oo) = — 3^/Ai so that A vanishes here also. 
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FIG. 4: (Color online) Solution plot for case Ai = for different initial conditions. 

VI. LINEARIZING AND CANONICAL TRANSFORMATIONS 

At this stage one may also try to trace the reason for the existence of the above type 
of remarkable oscillatory behaviour for Eq. (j2J). Under the transformation x = jp, Eq. (J2J) 
transforms to a linear third order differential equation, 



from which also one can trace the solution (|13|). Further, it gets transformed into a free 
particle equation X" = 0, where prime denotes differentiation with respect to a new variable 
r, with the point transformation 



w + Xiw = 



(25) 




3Al y/ — \ x x 




3x/^A7 x 



(26) 
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Even then it is very unexpected and unusual to realize such nonlinear oscillator systems 
having simple properties similar to that of linear oscillators. 

In addition to the above, one can also linearize the equation (J2J) through a nonlocal 
transformation, 

U = xe-zS x{t ' )dt ' . (27) 

Under this transformation Eq. (j2J gets modifieed to the form of a linear harmonic oscillator 
equation, 

i) + \ 1 U = 0. (28) 

Note that the above transformation (J27|) is valid for all the three cases, namely, Ai > 0, Ai = 
and Ai < 0. For Ai > 0, obviously the solution of (|2*5jl is 



U = Asw(ut + 6), (29) 



where A and 8 are arbitrary constants and frequency u> = yAi, which is independent of 
the amplitude. Consequently, from the forms of the general solutions (}2"9*|) and (JT3J) we can 
identify the obvious canonical transformation 

* = TZTEp' (30) 

1 3Ai r 

where P = U, so that using eq. (|15j) for the canonically conjugate momentum p, we have 

p = P(l-g^P). (31) 

It is straightforward to check that when U and P are canonical so do x and p (and vice versa) 
and that the Hamiltonian H in eq. (j!6|) can be rewritten as the standard linear harmonic 
oscillator Hamiltonian 

H= -{P 2 + \iU 2 ). (32) 

More interestingly, even for Ai < 0, making use of similar argument for the aperiodic 
bounded solution (jl9j) . we can write the canonical transformation 

x = ^— , p = P(l + -Lp) (33) 
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where P = U, so that the Hamiltonian (fTtjj) (for Ai < 0) is mapped on to the unbounded 
(wrong sign) linear harmonic oscillator 

H= 1 -{P 2 -uj 2 U 2 ), u=y^\. (34) 

Finally for the Ai = case, we can identify the canonical transformation 

3P kU 2 
X = W P= 6~' (35) 

so as to transform to a freely falling particle Hamiltonian H = — \J\U. However, we 
suspect that there may exist a canonical transformation for this specific case Ai = 0, which 
will take the corresponding form of Eq. (J2J) to a free particle, though we could not succeed 
to find it so far. 

Of course, it is well known that linearization in general does not ensure preservation 
of properties of linear systems in the nonlinear case. The classical example is the Cole- 
Hopf transformation through which the nonlinear Burgers equation is transformed into the 



linear heat equation ^j]. This is also true in the case of various soliton equations such as 
Korteweg-de Vries, sine-Gorden, nonlinear Schrodinger and other equations which are all 
linearizable in certain sense. Nevertheless in the present example of Eq. the amplitude 
independence of frequency does indeed get preserved. These considerations also raise the 
question of identifying the most general class of nonlinear dynamical systems of the form 
x = f(x,x,t) that has solutions whose frequency is independent of the amplitude as in the 
case of the linear harmonic oscillator, a point which is under study presently. 

VII. GENERALIZATION AND OTHER EXAMPLES 
A. Generalization 

We also observed that the above exceptional properties admitted by Eq. (J2J) is also 
common for a class of nonlinear systems. In particular we find that the modified Emden 
type equation with linear term and constant external forcing 

x + (ktx + k 2 )x + ^-x 3 + ^-^x 2 + \ lX + A 2 = 0, (36) 
9 3 

where k\,k2,Xi and A 2 are parameters, for which also one can obtain similar propertes as 
discussed above for Eq. (J2J) when the parameters satisfy the specific condition Ai = ~g- + ^f-- 
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To be specific, for the case -j^ — |- > 0, the explicit periodic solution takes the form 

A sin (ut + 5) 

xit) = 

l-(^)Acos(ut + Sy 



0<A<^, q=Jx 1 -% (37) 
where x(t) = x(t) + For ^j^- — ^ < 0, the system becomes aperiodic. 

B. Comparison with other nonlinear oscillators 

Finally, it is of interest to compare the dynamics of Eq. (j2J) with another interesting 
nonlinear oscillator of the form 

Xxx 2 lo 2 x . . 

^TTa^ + TW = °- (38) 

which possesses exact periodic solution x(t) = Asin(Qt + 5), where the frequency Q = 
Vi+aa 2 ' exn ibiting the characteristic amplitude-dependent frequency of nonlinear oscilla- 
tors inspite of the sinusoidal nature of the solution Eq. (J38j) . Equation (J38)) has a natural 
gen e r a llZation in three dimensions Q and these systems can be also quantized exhibiting 
many interesting features and can be interpreted as an oscillator constrained to move on 
a three- spher e llq . Several generalizations of these systems to A-degrees of freedom are 
also possible 17[. But all these nice properties are subject to the fact that the frequency is 
amplitude- dependent. It is only Eq. (j2J which admits the amplitude independence property 
of frequency for a nonlinear oscillator. Studies of such systems will have important impli- 
cations in developing nonlinear systems whose frequency remain unchanged even with the 
addition of suitable nonlinearity. Also, Eq. (J3"8j) when written as a system of two first order 
equations of the form (|18p. with 

Mx,y)=y and /a(s>y ) = _ 

the divergence of the flow function A = 7^ 0. On the other hand Eq. (J3%|) is also a 

Hamiltonian system Q| with the Hamiltonian given by 

H= l -[p\l + \x 2 )+ Ulx 1 



eventhough A ^ 0. Once again it is only 

A = - / r-^dt = 0. 

TJ l + Xx 2 
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This again confirms that for a conservative system it is necessary that the average of the 
divergence of the flow function vanish. 

VIII. CONCLUSION 

In this work, we have identified a nonlinear oscillator which exhibits unusual oscilla- 
tory properties: Purely harmonic oscillation whose frequency is completely independent of 
the amplitude and is the same as that of the linear harmonic oscillator. The conservative 
Hamiltonian nature of the system has been established, which also necessitates the gen- 
eralization of the definition of conservative systems in terms of the divergence of the flow 
function. Further, even for an apparently dissipative type/aperiodic (front like) type sys- 
tem, the Hamiltonian structure is preserved and the definition of the conservative system in 
terms of the divergence of flow function needs to be generalized. We believe identification 
of such nonlinear systems having the basic property of linear systems will have considerable 
practical application as the effect of higher harmonics is completely suppressed. Also it is 
of interest to consider the quantum mechanical version of the system and its higher 
dimensional generalizations as well as the effect of external forcing and additional damping, 
which are being pursued. 

We have also shown that Eq. (J2J) admits certain interesting geometrical properties. For 
example, under the invertible point transformation Eq. (J2J) gets transformed to the 

free particle equation whereas with the appropriate choice of canonical transformations one 
can transform it to a simple harmonic oscillator (or to an unbounded oscillator or to a 
freely falling particle depending on the value of the parameter Ai) equation. Interestingly, 
the very same nonlinear oscillator can also be transformed to a third order linear equation 
through a nonlocal transformation, while yet another nonlocal transformation transforms 
it to a linear harmonic oscillator. As far as our knowledge goes no such single nonlinear 
oscillator possesses such large class of transformation properties. We believe that exploring 
such nonlinear equations will be highly rewarding in understanding nonlinear systems and 
further work is in progress in this direction. 
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